Diffractive structure functions in nuclei 
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We calculate proton and nuclear diffractive structure functions in the IPsat (Kowalski- 
Teaney) dipole model. This parametrization has previously been shown to provide good 
agreement with inclusive F2 measurements and exclusive vector meson measurements at 
HERA. We discuss how the impact parameter dependence crucially affects our analysis, 
in particular for small 13. 

1 Introduction 

The large fraction of diffractive events observed at HERA shows that modern cohiders are 
approaching the nonhncar regime of QCD, where ghion saturation and unitarization effects 
become important. In Deep Inelastic Scattering (DIS) on nuclei the nonlinear effects are 
enhanced by the possibihty of interacting coherently with several nucleons simultaneously [5] . 
There are plans for several facilities capable of high energy nuclear DIS experiments, as the 
EIC [3] and LHeC 0] colliders. Due to the difficulty in measuring an intact recoil nucleus 
deflected by a small angle, diffractive eA collisions present an experimental challenge. But 
if they are successful, nuclear diffractive DIS (DDIS) would provide a good test of our 
understanding of high energy QCD. 

In the high energy limit DIS is best understood in the dipole frame, where the incoming 
virtual photon fluctuates into a quark- antiquark pair which then interacts with the target. 
The scattering amplitude is related to the correlator of two Wilson lines in the wavcfunction 
of the nucleus. In contrast to the language of coUinearly factorized parton distribution 
functions, in this formalism both inclusive and diffractive observablcs can be calculated 
from the same universal dipole cross section. This enables one to naturally use the dipole 
cross sections fitted to one process to predict observablcs in another one. In this talk we 
will review the results of our recent work [S] to apply this ideology to computing nuclear 
diffractive structure functions. Our emphasis is not on the most recent developments of 
high energy evolution equations, but the effects of a more realistic and consistent impact 
parameter dependence. This will lead us to discuss, in addition to the consequences of 
nuclear geometry on diffractive observablcs, the importance of the proton impact parameter 
profile used in the calculations. In this paper we shall first describe the dipole cross sections 
and calculation methods used and then summarize our results for nuclear DDIS. 

2 Method 

We decompose the diffractive structure function into different components in the standard 
way as 
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For the lowest Fock state of the virtual photon wavefunction, the qq dipole {F^qq and 
-^Lqq ) follow the treatment of [6j. At small /? (large mass of the diffractive system) 
the dominant contribution comes from higher Fock states. In this work we are interested 
in the finite experimentally relevant range of y/s and will only include the leading (in as) 
ggij-component of these. In different works this component has been evaluated in different 
limits, we shall here use the approach of [7j and interpolate between the small /3, large Nc 
formula used in [8] and the finite /9, large form used in [6]. We refer the reader to the 
references above for the detailed formulae. 

We use the "IPsat" dipole cross section parametrization introduced in Ref. [9] and ex- 
tensively studied in Ref. [10]. This is an impact parameter dependent dipole cross section 
that combines unitarization and the correct behavior of structure functions the logarith- 
mic large (i.e. small dipole size r) behavior of structure functions are achieved This is 
achieved using an eikonalized DGLAP-evolved gluon distribution function [11]. The dipole 
cross section is given by 
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where F is proportional to the DGLAP evolved gluon distribution 
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with both the coupling and the gluon distribution xg{x, Q^) evaluated at the scale /i^(r^) 

Several works on the subject (e.g. [T^l 
[6l [131 [H]) assume, explicitly or imphcitly, 
a factorizable 5^ dependence of the dipole 
cross section. 
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P = 0.9 
= 0.5 
P = (GBW)- 
P = (MS) 



which leads to an exactly exponential t- 
dependence of diffractive cross sections. 
The conceptual problem with the form 
Eq. (HJ) is that it cannot be a solution of the 
BK equation (unless the profile is a 6* func- 
tion, which would contradict the experimen- 
tally observed t-distribution). A factorized 
Gaussian profile for the proton dipole cross 
section, for example, does not approach the 
correct unitarity limit for b 0. This is the 

main motivation for including the impact parameter dependence in the saturation scale 
model (121); not as a factorizable prcfactor of the dipole cross section. For a proton the 
impact parameter profile in Eq. ^ is taken as T^bx) = Tp^br) ~ ^^'^ ^'^^ ^ 

nucleus T{bT) = X]iLi^p(^T — bxi), where the nucleon coordinates bxi are taken from a 
standard Woods-Saxon distribution [15]. The concrete consequence of this impact parame- 
ter dependence is that, in contrast to the factorized ansatz the different components of 



Figure 1: Contributions of different impact 
parameters to the inclusive and diffractive 
structure functions in the proton in the IPsat 
model. 
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Figure 2: Left: dependence on (3 of nuclear effects on the individual components of the 
diffractive structure function in a gold nucleus. Right: total /3-dependence of nuclear effects 
on the diffractive structure function. In both plots = 5 GeV^ and Xv = 10~^. 



the diffractive structure have different 6-dcpcndcnccs from each other and from the inclusive 
cross section. The gg-component is enhanced at small b (closer to the black disk limit), 
whereas the qqg-pait is dominated by larger b (because it vanishes in the black disk limit). 
This structure is illustrated in Fig. [T] 



3 Results 
3.1 HERA 

We compare our calculation to the HERA results on diffractive structure functions, measured 
using both using the rapidity gap method (ZEUS FPC [16] and HI LRG [17]) and by 
measuring the recoil proton (ZEUS LPS [18] and HI EPS [l^). Because the FPC and 
LRG data include events in which the proton has broken up, the cross-sections measured 
for the process ep —f eXY are larger than the one measured for the process ep — s- eXp. We 
scale down this data by a constant factor to correct for the proton dissociation contribution; 
the ZEUS FPC data by a factor of 1.45 and the HI LRG data by 1.23. These factors are 
different due to the different cuts on My, the mass of the proton dissociation system. For 
the combined datasct from ZEUS and HI data both with and without identified protons we 
get /Ndof = 1.3, with a coefficient as = 0.14 of the qqg term. This is the values of as that 
we shall use to evaluate nuclear diffractive structure functions in the next section. For the 
IPsat model the largest contribution to the comes from the rapidity gap method data 
at large (3. The fit to only the LPS (x^ = 0.5 IPsat) and EPS (x^ = 0.8) is much better. 
Considering just the LPS also accommodates a larger value of as = 0.21 with still < 1. 

The fit to HERA data is better with a smaller as than in Ref. [7] . Given the 6-dependence 
described previously this is to be expected. The factorized ^-dependence used in earlier 
calculations of the diffractive structure function such as Refs. [6l|7] forces the ggg-component 
to have the same impact parameter dependence as the gg-component. As discussed above, 
the qqg component is sensitive to larger impact parameters and is thus larger; in order to 
fit the same data this must be compensated by multiplying it with a smaller factor of as. 
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Figure 3: Left: dependence of the gold diffractive structure function on Q^. Right: depen- 
dence of the nuclear diffractive structure function on the mass number A for = 5 GeV^. 
In both plots xp = lO^-^. 



3.2 Predictions for nuclei 

In Fig. [2] we show the ratios of different components of the gold diffractive structure function 
to the proton one as a function of /3. The gg-components of the i^j^ are enhanced compared 
to A times the proton diffractive structure functions. This is to be expected, because of 
the fact that in a gold nucleus the dipole cross section is, on average over the transverse 
area, closer to the unitarity limit than the proton - it is "blacker" . The elastic scattering 
probability of a qq dipole is maximal in the "black disk" limit and the approach to it is quicker 
in a large nucleus. The qqg component, on the other hand, is suppressed for nuclei compared 
to the proton. This is due to the fact that in a nucleus the scattering amplitude is closer to 
the unitarity limit, when the qqg component vanishes. This leads to a nuclear suppression of 
the diffractive structure function in the small /3 region, where the qqg component dominates. 
The net result of the different contributions is that Fj^, for a large range in /3, is close to 
AF^^. In Fig. [21 we plot the total ratio as a function of /3 for different nuclei in the "non 
breakup" case. As expected from our prior discussion, one sees a strong enhancement with 
A for larger /3 and likewise, a stronger suppression with A at very small values of /3. 

In our formalism, if one requires that the nucleus stays completely intact, the average 
over the nucleon positions 6ti must be performed at the amplitude level. This is the case 
of coherent diffraction. If the nucleus is allowed to break up into color neutral constituents 
(referred to as incoherent diffraction), the average is performed at the level of the cross 
section. Measuring the intact recoil nucleus at such a small t experimentally at a future 
electron ion collider is challenging, so it is useful to consider both cases. A comparison of 
the "breakup" versus "non breakup" cross-sections can be seen in the left panel of Fig. [3] for 
the ratio of diffractive cross-sections as a function of . The results in Fig. [3] for the ratio 
of diffractive structure functions indicate that the diffractive cross-section in nuclei decrease 
more slowly for large than in the proton. This can be understood as a consequence of Qg 
being larger for nuclei and diffraction being much more sensitive to /Q^ than inclusive 
DIS. In the right panel of Fig. [31 the nuclear size A dependence of the longitudinal and 
transverse components of the diffractive structure function is shown for the "breakup" and 
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"non breakup" cases. In the "breakup" case, one sees a very weak A dependence. In the 
coherent "non breakup" case, one first notes that the diffractive structure function first 
decreases up to atomic numbers A ~ 10, before beginning to rise. As noted in Ref. [2], this 
is due to the typical scattering amplitude for small nuclei actually being smaller than for 
a proton because of the diluteness of the nucleus. This leads to a suppression of coherent 
diffraction. The "breakup" case, on the other hand, can only be enhanced in nuclei. For 
gold nuclei, the cross sections in the "non breakup" case are about 15% lower than in the 
"breakup" case. 

Because of the different nuclear modifications in inclusive and diffractive scattering, the 
fraction of diffractive events in an experiment depends on the detailed kinematics and exper- 
imental coverage. For moderate values of and large nuclei we expect a nuclear shadowing 
of the inclusive structure function by a factor 0.8 [5]. A typical nuclear enhancement of 
diffraction (at moderate values oi (3 > 0.2) is a factor of ~ 1.2. Combining these we expect 
CTD/ctot to be increased by a factor of 1.2/0.8 = 1.5 compared to the proton. Thus from a 
typical ep fraction of 15% we expect CTo/crtot to go up to 20% - 25% at an eA collider. 
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